Tail asymptotics for the total progeny of the 
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Summary. We consider a branching random walk on R with a killing 
barrier at zero. At criticality, the process becomes eventually extinct, and 
the total progeny Z is therefore finite. We show that P(Z > n) is of order 
(nln 2 (n)) _1 , which confirms the prediction of Addario-Berry and Broutin 
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1 Introduction 

We look at the branching random walk on M + killed below zero. Let b > 2 be a determinist 
integer which represents the number of children of the branching random walk, and x > be 
the position of the (unique) ancestor. We introduce the rooted b-axy tree T, and we attach 
at every vertex u except the root an independent random variable X u picked from a common 
distribution (we denote by X a generic random variable having this distribution). We define 
the position of the vertex u by 

S(u) := x + X v 

where v < u means that the vertex v is an ancestor of u. We say that a vertex (or particle) 
u is alive if S(v) > for any ancestor v of u including itself. 



The process can be seen in the following way. At every time n, the living particles split 
into b children. These children make independent and identically distributed steps. The 
children which enter the negative half-line are immediately killed and have no descendance. 
We are interested in the behaviour of the surviving population. At criticality (see below for 
the definition), the population ultimately dies out. We define the total progeny Z of the 
killed branching random walk by 

Z ■= #{u e T : S(v) >0Vv< u} . 
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Aldous [2] conjectured that in the critical case, E[Z\ < oo and E[Z\n(Z)] = oo. In pQ, 
Addario-Berry and Broutin proved that conjecture (in a more general setting where the num- 
ber of children may be random). As stated there, this is a strong hint that P(Z = n) behaves 
asymptotically like l/(n 2 ln 2 (n)), which is a typical behaviour of critical killed branching ran- 
dom walks. Here, we look at the tail distribution P(Z > n). We mention that the Branching 
Brownian Motion, which can be seen as a continuous analogue of our model, already drew 
some interest. Kesten [6] and Harris and Harris j5] studied the extinction time of the pop- 
ulation, whereas Berestycki et al. [3] showed a scaling limit of the process near criticality. 
Maillard [7j investigated the tail distribution of Z, and proved that P[Z = n) ~ n2 ^ - as 
expected. 

Before stating our result, we introduce the Laplace transform := £?[e tx ] and we 
suppose that 

• (p(t) reaches its infimum at a point t = p > which belongs to the interior of {t : 

4>{t) < oo}, 

• The distribution of X is non-lattice. 

The second assumption is for convenience in the proof, but the theorem remains true in the 
lattice case. The probability that the population lives forever is zero or positive depending 
on whether i?[e pX ] is less or greater than the critical value 1/6. In the present work, we 
consider the critical branching random walk which corresponds to the case i?[e pX ] = 1/6. 
For x > 0, we call P x the distribution of the killed branching random walk starting from x. 

Theorem 1.1 There exist two positive constants C\ and C 2 such that for any x > 0, we 
have for n large enough 



Hence, the tail distribution has the expected order. Nevertheless, the question to find an 
equivalent to P[Z = n) is still open. As observed in [I], in order to have a big population, a 
particle of the branching random walk needs to go far to the right, so that its descendance 
will be greater than n with probability large enough (roughly a positive constant). The 
theorem then comes from the study of the tail distribution of the maximum of the killed 
branching random walk. By looking at the branching random walk with two killing barriers, 
we are able to improve the estimates already given in pQ. 

The paper is organised as follows. Section 2 gives some elementary results for one- 
dimensional random walks on an interval. Section 3 gives estimates on the first and second 
moments of the killed branching random walk, while Section 4 contains the asymptotics on 
the tail distribution of the maximal position reached by the branching random walk before 
its extinction. Finally, Theorem 11.11 is proved in Section 5. 



(1 + x)e px 
nln 2 (n) 



< P X (Z >n)<C 2 



(1 + x)e px 
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2 



2 Results for one-dimensional random walks 

Let R n = Ro + Yi + . . . + Y n be a one-dimensional random walk and P x be the distribution 
of the random walk starting from x. For any k G R, we define (resp. r^T) as the first time 
the walk hits the domain (k, +00) (resp. (—00, k)), 

t£ := inf{n >0:R n >k}, 
r k := inf{n > : R n < k} . 

We assume 

(H) E[Yi] =0, 39, rj > such that E[e~ {e+ ^} < 00, E[e (1+r > )Yl ] < 00. 

All the results of this section are stated under condition (H). The results remain naturally 
true after renormalization as long as £?[e tyi ] is finite on a neighborhood of zero (and E\Yi] = 
0). Throughout the paper, the variables Cx, C2, • • • represent positive constants. We first 
look at the moments of the overshoot Uk and undershoot Lf. defined respectively by 

U k := S T + -k, 

k 

L k := k - S- . 

Lemma 2.1 There exists C 3 > such that E°[e Uk ] G [C 3 , 1 /C 3 ) /or any /c > and E°[e eLk ] G 
[C 3 ,l/C 3 }foranyk<0. 

Proof. This is a consequence of Proposition 4.2 in Chang [3]. □ 

The following lemma concerns the well-known hitting probabilities of R. 
Lemma 2.2 For any x > 0, 

(2.1) P x (r+ < To) = + o(l/k) . 

as k 00. Moreover, there exist two positive constants C4 and C5 stzcn inai, /or any real 
k > and any z G [0, k], we have 

(2.2) 04 £±i<p. W<T -)< ft £±i. 

Proof. Let > and a; G [0,k]. By Lemma |2.1[ we are allowed to use the stopping time 
theorem on (R n , n < min^o", t£)), and we get 

x = E X [R T +, r+ < r ] + E x [R r - , r Q < r+] . 

We can write it 

x = kP x {r+ <r ) + A 1 -A 2 
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where A\ and A 2 are nonnegative and defined by A\ := E x [Uk, < r ] and A 2 := 
E X [L , Tq < Tft]. Equivalently, 

(2-3) P X (r+<r ) = X - A l + A \ 

By Cauchy-Schwartz inequality and Lemma \2.1\ we observe that 

(A,) 2 < E*[U 2 k }P x {4 < To) < C,P x {rt < r ) . 
Since P x {t£ < r ) goes to zero when k tends to infinity, we deduce that 

lim Ax = 0. 

k— >oo 

By dominated convergence, we have also 

lim A 2 = E X [L ] 

k— >oo 

and E X [L ] < Cj by Lemma I2TT1 This leads to equation ( 12.1 ft since E[—S T -\ = x + E X [L ] 
Furthermore, we have < A\ < ^/~C§ and < A 2 < C-j. Therefore ( 12 .3p implies that 

px ( + . -\ < X + x + 1 

P (r k < r ) < — ^ < C 8 — . 

Similarly, 

x — \fCl 



P X (r: <r )> h 



We notice also that P x {r k < r ) > -P°(r fe h < r ). By (12. ip . there exists a constant C 9 > 
such that P°(r+ < r ~) > We get 

if x < \[Cl + 1 



P x,+ <T -\ > ! k+l ifx<yCg- 
< T ° J ^ \ C 10 ^ otherwise 



with C10 := ^: +2 - ^ implies that 



p-W < v) > o„f±l. 

Thus equation (12.21) holds with C 5 := C 8 and C 4 := Cn. □ 

Throughout the paper, we will write Afc(l) for any function such that 

< D x < A fc (l) < D 2 

for some constants D\ and D2 and large enough. The following lemma provides us with 
estimates used to compute the moments of the branching random walk in Sections 3 and 4. 
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Lemma 2.3 We have for any x > 0, 



(2.4) 
(2.5) 
(2.6) 



E 



E 



E 



k—x 



e^Xy*(J* + l),T+<7b 
e^e-^ + l),T+<T 

1=0 



A fc (l) 



A fc fi: 



A;' 

1 +x 



fc 2 ' 

A fc (l)(l + x) 



Proof. First let us explain how we can find intuitively these estimates. The terms of the 
sum within the expectation is big when Re is close to 0, and the time that the random walk 
spends in the neighborhood of before hitting level is roughly a constant. Moreover, by 
Lemma [2.11 we know that the overshoot and the undershoot L behave like a constant. 
From here, we can deduce the different estimates. In (I2.4p . the optimal path makes the 
particle stay a constant time near zero then hit level k which is of cost 1/k. In (12.51) . the 
particle first goes close to 0, which gives a term in (1 + x)/k, then go back to level k which 
gives a term in 1/k. Finally looking at (12.61) . we see that the particle goes directly to 0, 
which brings a term of order k because of the sum, and a term of order (l + x)/k because of 
the cost to hit before k. The proofs of the three equations being rather similar, we restrain 
our attention on the proof of (12. 4p for sake of concision. 

We introduce the function g(z) := e _z (l + z) and we observe that g is decreasing. Let also 



A:=E C 



e^£e-*(Jfc + l),r+<7b 



Let a > be such that P(Y 1 > a) > and P(Y 1 < —a) > 0. For ease of notation we suppose 
that we can take a — 1. For any integer i such that < i < k, we denote by the interval 
[i, i + 1), and we define 

Tf := M{n > 1 : R n G /<} , 
N{i) := #{n < min{r+,r -} : R n G /«} 

which respectively stand for the first time the walk enters Jj and the number of visits to the 
interval before hitting level k or level 0. We observe that 



A< 9(i + l)E°[e u "N(i),r+ <r ]. 



0<i<k 
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Let i be an integer between 1 and k — 1, and let z G [i, i + 1). We have 

P 2 (T, > min(r - r+)) > P 2 (P, < i2 x , W G [l,r -], R x < i) . 

We use the Markov property to get 

P\T t > min(r -,r+)) > E z [P h {r+ > r ) h=Rl , R, < i] . 

By Lemma I2~2l equation (12. ip (applied to —R), there exists a positive constant Cu such that 
ph ( T h > r )>C 12 /(l + h). This yields 

(2.7) P 2 (T. 4 > min(r -,r+)) > j^P(Ri < "1) =: Cjar^- 

When i < fc/2, (and 2 G [*,* + !)), we notice that 



< C u P z (t+ 2 <t ) 

< r 1 + 1 

where the last two inequalities come from Lemmas 12. II and I2.2L For i > k/2, we simply write 

^[e^,r+<r -]<sup^[e^]. 

fc>0 

Therefore, we have for any % < k, 

(2.8) E*[e v \ r+ < r ] < C M i±l. 

We obtain that for any integer i between 1 and k — 1, and any z G Ii, 
E* [e u *N(i), r+ < t»] 

< J^i 1 + n ) ( SU P pZ ( T i < min ( r fc + ' r o~)) ) SU P E * b Uk ' T k < min(r -, T t 

= ( 1 - sup P 2 (T; < min(r+, 7^))) sup P 2 [e^, r+ < min(r ~, T, 

< c^{ l + ifc 16 l + l 

(2.9) < C17 



A- 

(2 + 1) 3 
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by (12. 7p and (12.81) . We have to deal with the extreme cases i = and % > k — 1. For zG/q, 
we see that P z (Ti > min^o", t£)) > P(Y 1 < — 1), which yields by the same reasoning as 
before 

W [e^(0), r+ < r "] <C 1S \. 
Similarly, (\_k\ is the biggest integer smaller than k), 

E*[e u *N([k\),T+ <t ]<C 19 . 

Therefore, (12 .9p still holds for any integer i e [0, k), as long as C17 is taken large enough. By 
the strong Markov property, we deduce that 



E° [e u *N(i), r+ < r ] < C 17 P\T t < r A r } 
This gives the following upper bound for A: 



i + iy 



(2.10) 
In particular 



0<i<k 



k 



A<C l7 \ ^(z + l) 3 5 (z + l) = C 20 ^ 



0<i<fc 



with C 2 o := C\ 7 J2i>o^ + l) 3 fi'(^ + !)• This proves the upper bound of (12 .4p . For the lower 
bound, we write (beware that [4 > 0), 



e u *^Te- R '(R e + l),T+<Tv 



> P»(t? < r ). 



We apply ((2HJ) to get the lower bound of (pOL □ 



3 Some moments of the killed branching random walk 

For any a > and any integer n, we call Z n (a) the number of particles who hit level a for 
the first time at time n, 

Z n (a) := = n : t ~(«) > n - 1, r^(w) = n} 

where for any a, t~(u) is the hitting time of (—00, a) of the particle u. We notice that 
particles in Z n (a) can be dead at time n, but their father at time n — 1 is necessarily alive. 
Let also 

Z(a) : = Z„,(a). 

n>0 
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Similarly, for any k > a > 0, and any integer n > 0, we introduce 

Z n (a,k) := #{|u| = n : r _ (u) > ra- 1, r+(u) > n, t~(u) = n} , 
Z(a,k) := ^2z n (a,k). 

n>0 

In words, Z n (a, k) stands for the number of particles who hit level a at time n and did not 
touch level k before. 

We denote by S n = X + X 1 + . . . + X n the random walk whose steps are distributed like 
X. We define the probability Q y as the probability which verifies for every n 



(3.1) 



dPy\x ,..,x n ' (f)(p) n 
Under Q y , the random walk S n is centered and starts at y. 

Proposition 3.1 We have for any x > 0, and any a > 0, 

e p(fc-o) 



(3.2) 
(3.3) 

Besides, if x > a > 0, 
(3.4) 



E k [Z(a,k)} = A fc (l)- 
E k [Z(a,k) 2 ] = A fe (l) 



k ' 

p(2k~2a) 



E*[Z(a,k) 2 ] = A fc (l)(l + x) 



^p(k+x— 2a) 



Proof. Let y be any real in [0, k] and let a G [0,y]. We observe that 

E«[Z n (a, k)} = b n py(r > n - 1, r+ > n, r" = n) . 

The change of measure yields that 

£*[Z n (a,A;)] = e^me-^,To>n-l,T+>n,T-=n] 

= e^E y Q [^ a ' s "\ r - > n - 1, r+ > n, r" = n] . 

Summing over n leads to 

(3.5) Ey[Z(a, k)] = e^E y Q [e^, r~ < r+] . 

Suppose that y > k/2. We observe that 
^[e^,r-<r+] 

p( a-S 



Tfc /V , S - <a 

T k/2 



E h Q [* pLa \h=s 



'k/2 



T k/2 < T A t> S r^ > a 



k/2 
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We know by Lemma [27T1 that sup £<0 EQ[e pLe ] < C 22 . We deduce that 
% [e pLa , r- < r+] < C 22 (e p ^ k ^ + P y (r k/2 < r+ ' 

We use Lemma 12.21 (applied to Rg = k — Si) to see that for k greater than some constant 
K(a) (whose value may change during the proof), 

1 + k—y 



% [e pL % r- < r+] < G 



23" 



For y < k/2, we see that 



E y Q [e pL % r" < r+] < E% [e^] < C 22 



We deduce the existence of a constant C24 such that for any < a < y < k and any 
k > K(a), we have E y Q [e pL \ r~ < r+] < C 24 1± j= M . It yields by (J33J) that 



(3.6) 



E y [Z(a,k)] < C 24 e p{y - a)1 k ' 



k 



Since L a > 0, we get E v Q [e pLa , r a < t£] > Q y (r a <r£). By Lemma 

1 + k — y 



Q\r~ < r+) > G 



25- 



k — a 



Therefore, using ( 13. 5p . we get that 



(3.7) 



E y [Z(a, k)} > C 25 e p(y - a) 



1 + k—y 
k 



Equations (I3.6P and (13. 7p give (13. 2p by taking y = k. We turn to the proof of (13. 3p and (13. 4p . 

E y [Z{a,k) 2 } = ^E y \Z[a,k)Z n [a,k)\ 



n>0 



(3.8) 



n>0 |u|=n 

We decompose Z(a, k) along the particle u to get 

n-l 

Z(a,k) = l+^2z u *(a,k) 



£=0 



where ug is the ancestor of u at time I and Z Ut (a, k) is the number of descendants v of at 
time n which are not descendants of ug + i and such that n = r~(t>) < t£(v). In particular, 

E[Z u *{a,k)] = (b- 1) (E s ^ [E Sl [Z{a,k)l S l e [a,fc]] +P s( ^ ) (^i < «)) 
= (6-1)^,(1)^) 



e^- a ) 1 + fc ~ gl ,^€[a,fc] 



A- 



+ P(Y1 < a - S'(uf)) 



Afc(1)e P(SK)-a)i + i^M 

k 



if k > K(a) and S(ui) > a. This decomposition leads to 



Ey [Z(a,k),n = r-(u)<r+( 

-pa " 

^W^rE^ t^C* - + 1), n = r"(«) < r+(t*)] . 



£=0 



Then equation (13.81) becomes 

p -pa n 

E y [Z(a,k) 2 ] = A fc (l)— ^^^[e'^-ft+l) 



n = r„ 



n>0 £=0 
e p(y-a) 



V~ E E E Q [^- 5n) (^ - ft + 1), n = r- < r+] 



n>0 £=0 



(3.9) 



A,(r 



oP{y-°) 



£=0 



where we used the change of measure from P y to Q y defined in (13. ip . Take y = k. It implies 
that 



E k [Z(a, k) 2 ] 

p p{2k-2a) 



e ^*^ e - p(fc - s ^(ib - ft + 1), T~ < 



1=0 



We apply equation (12 .4p of Lemma [2731 for the walk i?^ := pik — ft) to get (13.31) . If we take 
y = x, we obtain 



pp(x+k-2a) 

E*[Z(a,k) 2 ] = A fe (l) ^ 



e^£e-^-^(A;-ft + l),r-<^ 



£=0 



and we apply (12.51) of Lemma 12.31 to complete the proof of (13 .4p . □ 

4 Tail distribution of the maximum 

We are interested in large deviations of the maximum M of the branching random walk 
before its extinction 

M := sup{S(u) : ueT such that S(v) > \/ v < u} . 
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To this end, we introduce 



H n (k) : = = n : r+(w) = n,r (u) > n} 

H(k) := Y,H n (k). 



n>l 



The variable H (k) is the number of particles of the branching random walk on [0, k] with 
two killing barriers which were absorbed at level k. 

Proposition 4.1 We have 

(4.1) E*[H k ] = A fc (l)e^- fc )^ 



(4.2) E x [Hl\ = A k {l)e pix - k) 



k 

1 + x 
k 



It shows that H k is strongly concentrated. Our result on the maximal position states as 
follows. 

Corollary 4.2 The tail distribution of M verifies 

P X (M> fc) = A fe (l)(l + x)^— 

Proof. The corollary easily follows from the following inequalities 

P X (M > k) < E x [H k ] 

and 

P(M>k) = P(H k >l)>^tf. □ 

We turn to the proof of Proposition 14.11 Since it is really similar to the proof of Propo- 
sition 13.11 we feel free to skip some of the details. 



Proof of Proposition 4 ■ 1 We verify that 

(4.3) E x [H k ] = e^E x Q [ e - pu \ r+ < r 1 . 
Since U k > 0, we deduce that 

(4.4) E x [H k ] < e^ x - k) Q x (r k + < r ) = A k (l)e p{x - k) ^^- 
On the other hand, observe that 

E x Q \e- plJ \ r+ < r "] > e~" M Q x {U k < M, r+ < r ) 

11 



We see that 



Q x (U k > M, r+ < to) 
< Q x (S T+ < k, r+ < Tq) sup Q° (Ut > M) + Q x (U k/2 > k/2) 

1 + X 



< 



k 



-e(M) + o(l/k) 



for some e(M) which goes to zero when M goes to infinity by Lemma [2.11 Therefore 

1 + x 



(4.5) 



Q x (U h < M, r+ < r ") > C, 



26" 



for M large enough. Equations (j4.3j) . (I4.4p and fl4.5[) give (I4.ip . We look then at the second 
moment of Hk- As before (see (13. 9p ). we can write 



E x m\ = A fc (i; 



3 p(x-fe) 



A- 



£=0 



We apply (12. 6 p of Lemma 12.31 to complete the proof. □ 



5 Proof of Theorem 11.1 

Proof of Theorem ! 1. 11 lower bound. Let a G (0, x). We observe that 

P X {Z >n)> P X {M > k)P k {Z{k, a)>n). 

By Proposition 13. \\ there exists a constant u > such that E k [Z(k,a)] < fie pk /k when fc is 
large enough. Let k be such that fie pk /(2k) = n. Then k = Mn(n) + o(ln(n)), and we get 
by Corollary |42] 

By the choice of k, we notice that 



P X (M >k)> C 27 { t 2 



nln (n) 



P fe (Z(ifc, a) > n) > P fc f Z(Jfe, a) > 



E[Z(k,a)} 



Thus Paley-Zygmund inequality leads to 



P k (Z(k,a) >n)>- 



lE k [Z(k,a)f 



AE k [Z(k,a 



21 ' 



Prop osit ion 13.11 shows then that P k (Z(k,a) > n) > C 2 s > 0. Therefore, 

P X (Z > n) > C j±±^- 
nln In) 
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with C29 = C27C28/4, which proves the lower bound of the theorem. □ 



We turn to the proof of the upper bound. We recall that Z(0) represents the number of 
particles who hit the domain (— 00, 0). 

Proof of Theorem : upper bound. First, we notice that Z(0) = 1 + (b — 1)Z. Indeed, Z(0) 
is the number of leaves of a tree of size Z + Z(0), in which any vertex has either zero or b 
children. Therefore 

P X (Z >n)=P x (z(0) > 71 ~ ' 



6-1 

Hence it is equivalent to find an upper bound for P x (Z(0) > n). For any k, we have that 

P x {Z{0)>n) < P X [M < k, Z(0,k) > n) + P X {M > k) 
< P x (Z(0,k)>n)+P x (M>k). 

By Markov inequality, then Proposition 13. 1\ we have 

E x \Z^kf\ ^ n n , . px e" k 



P x (Z(0, k)>n)< 1 Y ' J < C 30 (l + x)e 



n 2 " ' k 3 n 2 



Therefore, by Corollary 14. 2\ we have for k large enough 



P x (Z(0)>n)<C 31 (l + xy x [ ¥ ^^ 



e pk e~ pk 



Take k such that e pk /k = n. We verify that 

e pk e - P k 1 1 



k 3 n 2 k p 2 nln 2 (n) 
which gives the desired upper bound. □ 



(1 + 0(1)) 
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